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Abstract 

An  investigation  of  scattering  from  surface  cracks  has  been  conducted. 

In  particular,  the  change  in  the  reflection  coefficient  of  a  Rayleigh  wave 
incident  on  a  surface  indentation  crack  has  been  measured  as  the  sample  is 
stressed  to  fracture.  The  acoustic  measurements  have  been  correlated  with  the 
stable  crack  extension  that  precedes  final  failure.  The  crack  extension 
behavior  in  as-indented  and  annealed  specimens  was  found  to  be  appreciably 
different.  The  cracks  in  the  annealed  samples  exhibited  partial  crack  closure 
but  little  stable  extension,  whereas  the  cracks  in  the  as-indented  samples 
displayed  both  crack  closure  and  irreversible  crack  growth.  This  behavior  has 
been  rationalized  by  invoking  concepts  based  upon  the  residual  stresses 


created  by  indentation. 
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I.  INTRODUCTION 


The  use  of  surface  acoustic  waves  to  detect  and  characterize  surface 
cracks  in  brittle  materials  has  been  described  in  previous  studies. The 
basic  technique  devised  in  these  studies  requires  the  excitation  of  a  Rayleigh 
wave  on  the  surface  of  the  brittle  solid  and  the  subsequent  observation  of  the 
acoustic  wave  reflections  from  the  surface  located  cracks.  The  first  studies* 
were  performed  on  glass  samples  containing  surface  cracks  introduced  by 
indentation  and  then  extended  to  a  larger  size  by  application  of  a  sub- 
critical  stress.  The  acoustic  measurements  were  conducted  with  a  small  stress 
applied  to  the  specimen  in  order  to  eliminate  crack  closure.  These  studies 
indicated  that  acoustic  measurements  performed  in  the  long  wavelength  regime 
provided  precise  predictions  of  the  failure  strength  of  the  glass  specimens. 
Subsequent  studies  were  conducted  on  S13N4  specimens.^  The  surface  cracks 
subject  to  acoustic  measurements  in  these  studies  were  created  by  direct 
indentation  without  introducing  an  intermediate  extension  procedure.  The 
final  fracture  stresses  obtained  in  the  test  specimens  could  be  adequately 
predicted  from  the  acoustic  measurements,  but  only  by  invoking  an  effective 
fracture  toughness  smaller  than  the  actual  toughness  of  the  material  by  ~30% 
(see  Appendix).  Additionally,  the  acoustic  estimates  of  tfTe  crack  radius  were 
invariably  smaller  than  the  subsequently  measured  crack  radii  (at  the  fracture 
criticality)  by  a  factor  of  ~2.5  .  Cracks  of  this  type  are  typical  of 
machining-induced  cracks,-*  and  hence,  an  adequate  comprehension  of  such 
discrepancies  is  needed  before  applying  surface  acoustic  wave  techniques  to 
the  prediction  of  failure  of  components  subject  to  surface  finishing  by 
machining. 
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Recent  studies  of  indentation  fracture^"'  have  demonstrated  that  the 
cracking  behavior  is  dominated  by  the  residual  field  created  during  the 
indentation  process.  This  residual  field  causes  indentation  cracks  to  extend 
when  an  external  stress  is  applied.®  The  initial  extension  is  stable; 
instability  commences  when  the  crack  attains  a  length  ~2.5  times  the  initial 
length.  This  residual  stress  effect  is  almost  certainly  the  origin  of  the 
discrepant  acoustic  predictions.  The  phenomenon  is  examined  in  the  present 
study  by  monitoring  the  acoustic  scattering  from  a  surface  crack  as  an 
external  stress  is  applied  to  the  test  specimen. 

The  residual  stress  can  be  essentially  eliminated  by  annealing  at 
temperatures  (~1200°C  for  Si 3N4)  at  which  dislocation  climb  and  annihilation 
can  proceed.  Correctly  annealed  specimens  should  thus  exhibit  no  stable  crack 
extension  and  an  enhanced  fracture  strength.  Hence,  the  influence  of 
annealing  on  both  the  acoustic  scattering  and  the  failure  prediction  is  also 
examined. 

Finally,  implications  for  predicting  the  failure  of  machined  components 
from  surface  acoustic  wave  measurements  are  discussed.  Emphasis  is  placed 
upon  the  choice  of  an  optimum  measurement  procedure  and  the  role  of  annealing. 

II.  ACOUSTIC  SCATTERING  FROM  SURFACE  CRACKS 

A  general  theory  of  scattering  from  flaws^*^  forms  the  basis  for  the 
present  study.  The  scattering  configuration  is  shown  in  Fig.  1.  The  quantity 
of  interest  is  the  relative  amplitude  of  the  acoustic  signal  scattered  by  the 
flaw.  This  quantity  is  given  by  the  reflection  coefficient  S21  and  is 
defined  as  the  ratio  of  the  amplitude  of  the  signal  reflected  from  the  flaw. 


(received  by  transducer  2)  to  the  amplitude  of  the  incident  signal,  Aj 
(transmitted  by  transducer  1),  at  the  terminals  of  the  transducers.  Expressed 
in  terms  of  the  input  power  to  the  transmitting  transducer,  PR  ,  the  Rayleigh 
wave  displacement  field  (when  the  receiving  transducer  is  used  as  the 
transmitter)  uj^'  ,  and  the  prior  stress  in  the  vicinity  of  the  flaw, 
o^1)  ,  the  reflection  coefficient  for  a  void  is  given  by 
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Here,  the  integral  is  taken  over  the  entire  surface  of  the  void,  Sc  .  For  a 
Rayleigh  wave  normally  incident  on  a  half-penny  shaped  surface  crack  of 
radius  a  located  in  the  x-y  plane  (Fig.  2),  Eq.  (1)  becomes 
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where  auz  is  the  discontinuity  in  displacement  across  the  crack  and  ozz  is 
the  prior  normal  stress.  The  integral  is  taken  over  just  the  semi-circular 
area,  S  .  The  theory  is  strictly  valid  only  in  the  low  frequency  regime 
(i.e.  where  the  maximum  depth  to  which  the  crack  extends  below  the  sample 
surface  is  required  to  be  much  less  than  an  acoustic  wavelength).  The 
scattering  from  surface  cracks  is  also  influenced  by  imaging  at  the  surface 
(which  increases  the  value  of  the  stress  intensity  factor)  and  the  variation 
with  depth  of  the  Rayleigh  wave  stress  field.  These  influences  are 
incorporated  into  the  analysis  by  adopting  the  results  of  Budiansky  and 
O'Connell^  which  permit  the  surface  integral  in  Eq.  (2)  to  be  expressed  in 
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terms  of  a  contour  integral 
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where  v  is  Poisson's  ratio,  E  is  Young’s  modulus,  Kj  is  the  mode  I 
stress  intensity  factor,  and  p(r)  is  defined  for  a  given  point  r  on  the 
crack  circumference  C  as  the  normal  to  the  tanrent  line  at  r  which  also 
passes  through  the  origin  of  the  crack  coordinate  system.  The  angular 
dependence  of  Kj  caused  by  the  surface  imaging  forces  may  be  approximated 
from  the  results  of  Smith,  Emery,  and  Kobayashi^  for  a  half-penny  shaped 
surface  crack  in  a  beam  of  thickness  2c  subject  to  a  bending  load.  The 
applied  stress  then  takes  the  linear  form 


°zZ(y)  =  -  y/c) 


(4) 


where  8  is  a  constant  and  y  is  the  distance  from  the  sample  surface 
(Fig.  2).  The  corresponding  stress  intensity  factor  is  given  by 


Kj (6)  =  2/a7?  8[to(0)  -  (a/c)  ♦j(e)]  (5) 

where  the  functions  <|>q(6)  and  ^(e)  have  been  numerically  evaluated 
(Fig.  3).  These  solutions  are  used  to  provide  an  approximate  surface  wave 
scattering  result  by  adopting  a  linear  approximation  to  the  Rayleigh  wave 
stress  field  and  deducing  effective  values  for  the  constants  8  and  c  .  The 
normalized  reflection  coefficient  wk|Sjj|  for  a  semi-circular  crack  obtained 


in  this  manner  is  plotted  as  a  function  of  the  normalized  crack  depth, 

2ira/A  ,  in  Fig.  4.  Here  k  =  2n/\  is  the  propagation  constant  of  the 
Rayleigh  wave,  and  w  is  the  width  of  the  acoustic  beam  at  the  location  of 
the  crack. 

III.  EXPERIMENTAL  RESULTS 

A  commercial  hot-pressed  silicon  nitride  ( NC- 132)  was  used  for  the 
present  study.  The  samples  were  in  the  form  of  7.6  cm  x  2.6  cm  x  0.64  cm 
plates  with  polished  surfaces.  A  crack  was  introduced  into  the  center  of  each 
plate  surface  (normal  to  the  longer  edge)  by  Knoop  indentation  at  50N 
load.  Acoustic  measurements  were  made  at  a  frequency  of  8.5  MHz  , 
corresponding  to  an  acoustic  wavelength  of  ~680  urn  .  Two  wedge  transducers 
in  the  configuration  shown  in  Fig.  1  (with  a  =  12.5°)  were  used  to  excite 
and  receive  the  acoustic  signals.  The  experimental  values  for  the  reflection 
coefficient  we^e  corrected  for  diffraction  loss,  the  surface  wave  conversion 
efficiencies  of  the  transducers,  and  for  the  inclination  of  the  transmitted 
and  reflected  acoustic  beams  to  the  crack  normal. 

Measurements  were  performed  on  both  annealed  and  as-indented  samples. 

One  sample  was  annealed  in  air  and  another  in  a  vacuum  for  an  annealing  period 
of  about  six  hours  at  a  temperature  of  1200°C  .  An  additional  three  samples 
were  tested  in  the  as-indented  state.  In  all  cases,  the  change  in  the 
acoustic  reflection  coefficient  was  monitored  while  the  samples  were  slowly 
stressed  to  fracture  in  a  flexural  mode.  Typical  variations  of  the  reflection 
coefficient  (expressed  in  terms  of  a  predicted  crack  radius  -  see  Section  IV) 
with  the  magnitude  of  the  applied  flexural  stress  a'-e  plotted  in  Fig.  5. 


-6- 


Substantially  different  behavior  is  apparent  in  the  as-indented  and  annealed 
specimens.  The  as-indented  samples  exhibit  a  large  change  in  the  reflection 
coefficient  as  the  flexural  stress  is  slowly  increased  to  failure.  Also,  the 
change  is  largely  irreversible,  as  manifest  in  a  permanent  increase  in 
following  the  load  removal  prior  to  failure  (Fig.  5a).  This  behavior  suggests 
stable  crack  extension  during  the  application  of  an  external  stress, 
qualitatively  consistent  with  the  residual  stress  notions  of  indentation 
fracture.®  The  existence  of  stable  crack  extension  was  substantiated  by 
observations  of  the  fracture  initiation  site  on  the  fracture  surface 
(Fig.  6a).  These  observations  indicate  an  inner  crack  front  location  that 
coincides  with  the  initial  radius  ag  of  the  indentation  crack  and  an  outer 
crack  front  with  a  radius  ac  »  2.5ag  .  The  outer  front  is  presumed  to  be 
the  location  of  the  crack  at  the  fracture  criticality,  as  suggested  by 
observations  on  other  materials. 

The  annealed  specimens  (Fig.  5b)  exhibited  appreciably  larger  final 
levels  of  the  fracture  stress  (the  elevation  of  the  fracture  stress  is  similar 
to  that  measured  in  prior  studies).^  The  vacuum-annealed  specimen  exhibited 
small  reversible  changes  in  the  acoustic  scattering  during  stressing,  a 
phenomenon  that  will  be  attributed  to  the  gradual,  reversible  separation  of 
the  crack  surfaces  (see  Section  IV. B).  Observations  of  the  fracture  origin 
(Fig.  6b)  revealed  closely  spaced  crack  fronts,  indicating  that  stable  crack 
extension  had  been  largely  eliminated  by  the  annealing  treatment.  Some  stable 
extension  may  have  occurred  just  prior  to  instability  in  the  range  of  rapidly 
varying  crack  length  indicated  in  Fig.  6b. 
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The  specimen  annealed  in  air  (Fig.  5c)  exhibited  an  additional  increment 
in  the  scattering  coefficient  at  about  one  third  of  the  fracture  stress.  This 


effect  may  be  attributed  to  the  fracture  of  the  oxide  film  that  forms  within 
the  crack  when  annealing  is  conducted  in  an  oxidizing  environment. 

IV.  DISCUSSION 
A.  As-Indented  Specimens 

The  acoustic  scattering  and  the  fracture  response  exhibited  by  as- 
indented  specimens  are  both  influenced  by  specific  characteristics  of  the 
indentation  plasticity  and  the  indentation  morphology.  The  important 
relationship  between  scattering,  fracture,  and  indentation  can  be  established 
by  invoking  the  appropriate  indentation  parameters.  The  acoustic  scattering 
from  an  indentation  crack  can  be  interpreted  by  noting  that  the  residual 
stress  field  (attributed  to  the  elastic  constraint  of  the  indentation 
plasticity)^  allows  the  crack  surfaces  in  the  elastic  region  to  remain  fully 
separated  when  the  indentation  load  is  removed.  However,  the  crack  terminates 
in  the  vicinity  of  the  elastic/plastic  boundary  (Fig.  7a).  The  crack  thus 
exhibits  the  morphology  of  a  semi -annul us,  with  the  inner  radius  dictated  by 
the  radius  of  the  plastic  zone.  This  crack  geometry  is  not  amenable  to 
acoustic  scattering  analysis  using  conventional  analytic  techniques.  An 
approximate  solution  is  thus  developed  for  preliminary  crack  size 
estimation.  The  solution  considers  the  open  region  of  the  crack  to  be 
elliptical  in  section,  with  the  major  axis,  2a  ,  given  by  the  crack  diameter 
and  the  minor  axis,  2b  ,  given  by  the  difference  between  the  crack  depth  and 
the  plastic  zone  depth  (Fig.  7b).  In  essence,  therefore,  the  crack  is 


considered  to  be  held  fully  closed  at  the  specimen  surface  by  virtue  of  the 
integrity  of  the  plastic  zone.  The  scattering  induced  by  an  elliptical  crack 
subject  to  a  linearly  varying  stress  (the  linearized  Rayleigh  wave  stress 
field)  should  thus  provide  the  approximate  solution  of  present  interest.  The 
stress  intensity  factor  that  characterizes  this  configuration  is  given  by^ 
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where  K(k)  and  E(k)  are  complete  elliptical  integrals  of  the  first  and 
second  kinds,  respectively,  with  k^  =  1  -  (b/a)^  and  (k')2  =  (b/a)^  . 

My(e)  and  ML(e)  represent  the  stress  intensity  magnification  factors 
(Fig.  8)  for  an  elliptical  crack  subject,  respectively,  to  a  uniform  applied 
stress  and  a  linearly  varying  applied  stress,  and  Cj  and  C2  are  constants 
that  characterize  the  linearized  Rayleigh  wave  stress 
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The  stress  intensity  factor  can  be  inserted  into  Eqs.  (2)  and  (3)  in  order  to 
compute  the  scattering  coefficient.  The  result  for  the  normalized  reflection 
coefficient,  wk|Snl  »  as  a  function  of  the  normalized  crack  depth, 

2ir(b  +  h ) /A  ,  is  given  by  the  solid  curve  in  Fig.  4.  A  comparison  of  the 
acoustic  estimates  for  the  initial  crack  depth,  Cq  ,  and  the  depth  at 
fracture,  Cm  ,  with  the  corresponding  values  measured  optically  is  given  in 
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Table  I.  From  the  acoustic  estimates  for  the  initial  crack  depth,  a 
prediction  for  the  fracture  stress  of  each  sample  may  further  be  made  by  using 
the  relation 


where 


Acoustic 

Or 


kI  surf 


(8) 


kI | surf  =  1*22kl|bulk  =  1 . 22 ( 2/a7" )  (9) 

and  Kj£  =  3.55  MPa/m  These  predictions  for  the  fracture  stresses  are 
compared  with  the  actual  fracture  stress  values  in  Table  II.  Considering  the 
approximate  nature  of  the  scattering  analysis,  both  comparisons  yield 
remarkably  good  agreement. 

The  fracture  process  is  strongly  influenced  by  the  magnitude  of  the 
residual  stress  induced  by  the  plastic  zone.  The  extension  of  an  indentation 
crack  under  the  combined  influence  of  an  external  stress  and  the  residual 
indentation  stress  has  been  predicted  for  axisymmetric  indentations  by 
invoking  residual  stress  concepts.6  The  direct  applicability  of  these 
predictions  to  elongated  indentations,  such  as  Knoop  indentation,  has  not  yet 
been  verified  in  detail.  However,  it  is  instructive  to  compare  the 
axisymmetric  predictions  with  the  acoustic  measurements  of  crack  extension. 

The  following  expression  relates  the  equilibrium  crack  radius,  a  ,  to  the 
external  stress  o 

oo 
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(10) 


where  Kjq  is  the  material's  fracture  toughness,  P  is  the  indentation  load, 
xr  is  a  parameter  (related  to  material  hardness  and  elastic  modulus)  which 
determines  the  intensity  of  the  residual  field,  and  n  is  a  material 
independent  crack  geometry  constant.  According  to  Eq.  (10),  the  crack  length 
increases  stably  with  applied  stress  from  the  initial  value  (given  by  Eq.  (10) 
with  a  =  0) 

a0  =  (xrP/KIC)2/3  (11) 

to  the  failure  point  (given  by  do^/da  =0  in  Eq.  (10)) 

'  <4xrP/KIC)2/3  (12) 

om  ■  (27/256)1/3  [KIC/(,Q)1/2][KIC/xrP]1/3  (13) 

The  parameters  Kj^(nn)^^  and  xrP/KjQ  required  for  plotting  Eq.  (10)  can 
be  conveniently  calibrated  from  the  calculated  (from  acoustic  measurements) 
initial  crack  length  and  the  measured  fracture  stress.  Setting  ag  =  96  ym 
in  Eq.  (11)  gives  XpP/Kjc  =  9.2  x  10"V^  ,  while  the  measured  strength 
om  =  240  MPa  in  conjunction  with  Eq.  (13)  gives  K^/fTiTT  =  4.95  MPa/m  .  The 
predicted  dependence  of  crack  extension  on  applied  stress  is  shown  in 
Fig.  5a.  The  prediction  correlates  well  with  the  acoustic  measurements  of 
changes  in  crack  length. 


B.  Annealed  Specimens 

Removal  of  residual  stress  by  annealing  results  in  the  development  of 
crack  surface  contacts  over  an  appreciable  fraction  of  the  surface  (Fig.  9). 
Complete  closure  over  the  contacting  region  is  prevented  by  "asperities"  along 
the  crack  surface.  The  asperities  comprise  either  grains  dislodged  by  crack 
bifurcation  or  exposed  grains  in  which  residual  stress  (due  to  thermal 
contraction  anisotropy)  has  been  relaxed  by  the  passage  of  the  crack,  thus 
causing  the  grains  to  displace.  Compressive  forces  must  develop  at  the 
contact  locations  due  to  elasticity  of  the  surrounding  material.^  The  forces 
are  maximum  close  to  the  crack  tip  and  decrease  to  zero  at  some  distance  a.,- 
(Fig.  9).  The  area  subject  to  zero  contact  forces  defines  a  non-contacting 
crack  radius  ai  .  The  compressive  forces  allow  the  contacts  to  be  retained 
in  the  presence  of  an  acoustic  stress.  Hence  the  scattering  from  the 
contacting  region  of  the  crack  surface  is  negligible*  and  the  acoustic 
scattering  amplitude  only  provides  a  measure  of  the  non-contacting  crack 
radius  (Fig.  9). 

The  compressive  tractions  at  the  asperities  are  relaxed  by  application  of 
an  external  tensile  stress.  An  analysis  of  the  crack  separation  process^ 
indicates  that  separation  initiates  from  the  periphery  of  the  initial  non¬ 
contacting  zone  and  extends  as  the  external  stress  increases.  Full  crack 


Scattering  occurs  from  the  regions  between  the  asperities  within  the  contact 
area.  However,  since  these  scattering  centers  act  as  small  cracks  (with  radii 
approximately  equal  to  several  grain  diameters)  and  the  scattering  amplitude 
is  (to  lowest  order)  proportional  to  a^  ,  the  total  scattered  amplitude  from 
the  contacting  area  is  small  compared  with  that  from  the  non-contacting  area. 
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surface  separation  occurs  immediately  prior  to  unstable  crack  extension.  The 
crack  separation  process  is  reversible,  and  this  process  fully  accounts  for 
the  changes  in  acoustic  crack  length  that  occur  during  a  breaking  test. 

V.  IMPLICATIONS  AND  CONCLUSIONS 

A  combined  acoustic,  fracture  mechanics  study  has  demonstrated  that  long 
wavelength,  surface  acoustic  wave  scattering  measurements  can  provide  good 
estimates  of  the  dimensions  of  surface  located  cracks  and  that  the  crack 
dimensions  can  be  clearly  related  to  the  ultimate  failure  strength.  The 
connection  between  the  failure  condition  and  the  initial  crack  size  estimate 
depends  upon  the  state  of  residual  crack  stress  induced  by  the  machining  (or 
indentation)  process.  The  most  severe  residual  stress  exists  immediately 
following  the  machining  increment  that  generated  the  crack  subject  to  acoustic 
investigation.  Subsequent  machining  passes  will  remove  portions  of  the 
plastic  zone  and  successively  reduce  the  residual  stress.  However,  the  state 
of  residual  stress  at  a  specific  machining-induced  crack,  identified  by 
acoustic  inspection,  will  generally  be  unknown.  A  failure  prediction 
procedure  that  recognizes  the  most  detrimental  residual  stress  condition  is 
thus  recommended. 

The  most  consistent  acoustic  crack  size  estimation  can  be  conducted  in 
the  presence  of  the  residual  stress  because  the  crack  surfaces  are  fully 
separated.  Therefore  the  acoustic  measurements  should  be  conducted  on  the  as- 
machined  surface.  The  failure  strength  prediction  should  be  based  upon  the 
crack  size  estimation  procedure  pertinent  to  indentation  cracks,  using  the 
effective  value  of  the  initial  stress  intensity  factor  (Appendix)  that  fully 
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accounts  for  the  (residual  stress-induced)  stable  extension  prior  to 
failure.  This  procedure  will  provide  a  conservative  prediction  of  failure  for 
many  of  the  cracks  detected  on  a  machined  surface  because  some  of  the  residual 
stress  will  have  been  relaxed  by  partial  plastic  zone  removal.  However,  the 
unknown  state  of  the  residual  stress  necessitates  the  conservative  choice. 
Finally,  it  is  recommended  that  the  components  be  annealed  prior  to  use  in 
order  to  relieve  the  remnant  residual  stresses.  The  stress  relief  will 
provide  an  appreciable  margin  of  safety  for  most  cracks. 

At  this  juncture,  it  is  not  considered  appropriate  to  anneal  the 
component  prior  to  acoustic  evaluation  because  the  variability  in  the  crack 
closure  is  likely  to  provide  imprecise  estimates  of  the  actual  crack  length. 
However,  further  investigation  of  crack  closure  effects  may  suggest  effective 
procedures  for  generating  consistent  acoustic  measurements.  Finally,  it  is 
noted  that  the  strength  of  annealed  specimens  deduced  from  acoustic 
measurement?  on  machined  surfaces  may  be  a  consistent  underestimate,  whereupon 
the  failure  strength  prediction  could  be  elevated  above  the  minimum  level 
presently  recommended.  The  origin  of  such  an  underestimate  would  reside  in 
the  relative  influences  of  crack  closure  and  of  the  residual  stress  level  on 
the  relation  between  the  acoustic  scattering  amplitude  and  the  residual 
strength.  Again,  studies  of  crack  closure  effects  will  permit  such 
possibilities  to  be  fully  exploited. 
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APPENDIX 


Relation  Between  Failure  Strength 
and  Initial  Size  of  Indentation  Cracks 


Previous  studies  have  demonstrated  that  in  order  to  relate  the  final 
fracture  stress  of  indentation  flaws  to  the  initial  flaw  size  using  the 
Griffith  strength  equation,  an  effective  toughness  smaller  than  the  actual 
material  toughness  must  be  invoked.  This  result  can  also  be  derived  from  a 
fracture  mechanics  analysis  which  takes  into  account  the  influence  of  residual 
stresses  induced  by  the  plastic  indentation  zone.  For  axisymmetric 
indentations  the  fracture  stress  can  be  written  (using  Eqs.  (11)  and  (13) 
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1/3 


'IC 


(irflan) 
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(A-l) 


Hence,  in  the  idealized  case,  the  effective  toughness  is  Kef.f  =  0.47  Kjq  . 

In  practice,  several  factors  can  increase  Keff  above  the  idealized  value. 
Firstly,  the  incidence  of  lateral  cracking,  or  chipping,  after  the  formation 
of  the  median  crack  diminishes  the  residual  field,  as  characterized  by  the 
parameter  xr  in  Eq.  (13).  This  effect  can  cause  Keff  to  exceed  the 
idealized  value.  Secondly,  in  materials  susceptible  to  subcritical  crack 
growth,  the  initial  flaw  size  can  increase  under  the  influence  of  the  residual 
stress,  during  the  time  period  between  the  indentation  and  strength  tests. 
Again,  the  consequence  is  an  increase  in  Keff  .  Finally,  the  influence  of 
the  elongated  geometry  of  the  Knoop  indentation  may  alter  the  numerical 
constant  in  Eq.  (A-l).  Typically,  for  Knoop  indentation  pre-cracks  in  Si 3N4  , 
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TABLE  I:  As-Indented  Samples 

Cg  =  initial  crack  depth 
Cm  =  crack  depth  at  fracture 


Acoustic 
Cg  (wro) 

Actual 

Cg  (urn) 

Deviation 

(Co) 

Acoustic 
Cm  ^m) 

Actual 

Cm 

Deviation 

97 

116 

16% 

184 

220 

16% 

94 

112 

16% 

191 

223 

14% 

96 

112 

14% 

185-191 

205 

7-10% 
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TABLE  II:  Fracture  Stresses  (As-Indented  Samples) 


Sample 

Acoustic 

Op  (MPa) 

Actual 
oF  (MPa) 

Deviation 

4 

262 

238 

10% 

5 

266 

239 

11% 

6 

263 

242 

9% 
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FIGURE  CAPTIONS 


A  schematic  of  the  geometry  considered  in  the  derivation  of  the 
reflection  coefficient,  S21  ,  for  a  flaw.  Here,  n  is  the  inward 
normal  to  the  flaw  surface,  Sc  ,  and  a  is  the  inclination  of  the  flaw 
to  the  transducers. 

Scattering  geometry  for  a  Rayleigh  wave  normally  incident  on  a  half-penny 
shaped  crack. 

Angular  variation  of  the  functions  ^g(9)  and  ^(0)  (after  Smith  et 
al.).11 

Comparison  of  the  results  for  the  normalized  reflection  coefficient, 
wk I  Si 1 1  ,  as  a  function  of  the  normalized  crack  depth  for  the  half-penny 
shaped  crack  theory  and  the  elliptical  crack  theory.  The  normalized 
crack  depth  is  defined  for  each  theory  by  the  quantities  2n a/ A  and 
2n (b  +  h)/x  ,  respectively. 

The  variation  of  the  crack  depth  (obtained  from  the  acoustic  scattering 
amplitude)  with  the  flexural  stress. 

(a)  As-indented  specimen:  note  the  irreversibility  of  the  crack  length 
indicated  by  the  permanent  change  following  unloading.  Also  shown  is  the 
crack  change  predicted  from  the  indentation  model; 

(b)  Vacuum-annealed  specimen  indicating  irreversibility; 

(c)  Oxygen-annealed  specimen. 

Optical  micrographs  of  fracture  origins: 

(a)  As-indented  specimens  indicating  an  initial  crack  front  much  smaller 
than  the  crack  front  at  instability; 


(b)  Vacuum-annealed  specimens  showing  a  small  zone  of  stable  crack 
extension. 
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7.  (a)  A  schematic  of  an  as-indented  crack  indicating  that  crack  closure  is 


exhibited  over  the  plastic  zone; 

(b)  A  subsurface  elliptical  crack  analog  of  (a)  used  for  acoustic 
scattering  analysis. 

8.  Stress  intensity  magnification  factors  M^.(e)  and  ML(e)  as  functions 
of  the  angle  o  for  the  elliptical  crack  configuraion  of  Fig.  8b, 
subject  to  a  uniform  applied  stress  and  a  linearly  varying  applied 
stress,  respectively. 

9.  The  crack  configuration  in  annealed  specimens  indicating  the  zone  of 
contact  and  the  compressive  forces  that  exist  at  asperities  in  the 
region. 
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ABSTRACT 


We  have  extended  our  earlier  investigation  of  scattering  from  surface  cracks.  In  particular,  we  have 
studied  the  change  in  the  reflection  coefficient  of  a  Rayleigh  wave  incident  on  a  half-penny  shaped  sur¬ 
face  crack  along  with  the  corresponding  change  in  the  acoustic  crack  size  estimates  as  the  cracked  sample 
is  stressed  to  fracture.  We  have  examined  in  this  manner  both  cracks  in  annealed  samples  and  as-indented 
cracks.  We  have  found  that  the  fracture  behavior  for  cracks  in  these  two  types  of  samples  differ  quite 
significantly,  with  the  cracks  in  the  annealed  samples  exhibiting  a  partial  crack  closure  characteristic 
and  the  cracks  in  the  as-indented  samples  displaying  both  crack  closure  and  crack  growth  effects. 


INTRODUCTION 


This  work  is  an  extension  of  our  earlier  study1 
aimed  at  the  establishment  of  procedures  for  locat¬ 
ing  and  characterizing  surface  cracks  in  structural 
ceramics.  The  basic  technique  we  have  been  using 
consists  of  launching  a  Rayleigh  wave  on  the  sur¬ 
face  of  the  ceramic  and  observing  the  reflections 
of  the  acoustic  surface  wave  from  the  crack.  The 
particular  type  of  crack  we  have  been  studying  is 
a  half-penny  shaped  surface  crack  introduced  at  a 
given  orientation  into  the  ceramic  surface  by  a 
Knoop  hardness  indentor.  In  a  previous  paper,1  we 
described  a  scattering  theory  valid  in  the  low 
frequency  regime  based  on  the  model  of  an  open  half¬ 
penny  shaped  crack.  This  theory  related  the  reflec¬ 
tion  coefficient  measured  in  the  experiment  with 
the  crack  size  and  fracture  stress  of  the  cracked 
sample.  The  theory  was  shown  there  to  give  predic¬ 
tions  for  the  fracture  stress  which  were  in  good 
agreement  with  experimental  results.  However, 
predictions  for  the  crack  size  were  observed  to  be 
considerably  less  accurate,  with  crack  size  esti¬ 
mates  for  unstressed  samples  being  smaller  by 
factors  of  two  to  three  than  the  actual  crack  sizes 
measured  after  fracture. 

Since  that  time  we  have  extended  the  theory 
and  correlated  our  results  with  what  would  be 
expected  from  fracture  mechanics  studies  of  cer¬ 
amics.  A  series  of  experiments  on  annealed  and 
unannealed  samples  has  demonstrated  the  power  of 
acoustic  techniques  for  studying  cracks  in  ceramics. 
The  work  indicates  that  cracks  in  annealed  samples 
tend  to  be  in  contact  over  most  of  their  cross- 
section  and  open  up  with  applied  stress.  At  the 
point  of  fracture,  a  very  slight  growth  in  crack 
radius  is  observed.  On  the  other  hand,  cracks  in 
unannealed  as-indented  samples  exhibit  partial 
closure  at  the  sample  surface.  Upon  application  of 
stress,  crack  growth  occurs,  with  crack  radii 
tending  to  increase  on  the  order  of  50%  .  Upon 
release  of  stress,  the  cracks  partially  close,  with 
their  effective  radii  decreasing  by  less  than  10%  . 
As  stress  is  further  applied,  further  crack  growth 


is  observed,  terminating  in  fracture  at  lower 
applied  stresses  than  for  equivalent  annealed 
samples.  This  hysteresis  effect  has  been  observed 
for  the  first  time  acoustically.  Both  the  initial 
crack  radius,  Cq  ,  and  final  radius,  Cm  ,  can  be 
measured  after  fracture.  The  results  obtained  can 
be  predicted  accurately  via  acoustic  techniques  and 
agree  well  with  the  predictions  of  earlier  theories 
and  experiments  by  Marshal  1 .2*3, 4 

An  important  result  of  our  work  is  the  demon¬ 
stration  that  it  is  highly  desirable  to  anneal 
ceramics  in  which  surface  cracks  are  likely  to  be 
present.  Annealing  inhibits  further  crack  growth 
by  relieving  the  residual  stresses  that  are  present. 
As  evidence,  we  note  that  when  unannealed  samples 
are  stressed,  the  crack  radius  increases  by  a 
factor  of  approximately  two.  Annealing,  however, 
causes  an  increase  in  the  fracture  stress  observed 
by  a  factor  of  approximately  1.5  . 

THEORETICAL  REVIEW 


A  general  theory  of  scattering  from  flaws 
developed  by  Kino  and  Auld^  forms  the  basis  for 
our  work.  The  scattering  configuration  considered 
is  shown  in  Fig.  1.  The  reflection  coefficient, 
S21  ,  is  defined  as  the  amplitude  ratio  of  the 
reflected  signal  from  the  flaw,  A2  ,  received  by 
transducer  2  ,  to  the  incident  signal,  A^  , 
transmitted  by  transducer  1  ,  at  the  terminals  of 
the  transducers.  In  terms  of  the  input  power  to 
the  transmitting  transducer,  P  ,  the  Rayleigh 
wave  displacement  field  when  the  receiving  trans¬ 
ducer  is  used  as  the  transmitter,  u V 2 )  ,  and  the 
applied  stress  in  the  vicinity  of  theJflaw  before 
the  flaw  is  introduced,  ,  S21  is  given  by 


;d 


for  the  case  where  the  flaw  is  a  void.  Here,  the 
integral  is  taken  over  the  entire  surface  of  the 
void,  S  .* 


Sc 

FLAW 


n  =  (  n  | , n2 , n3 ) 


Figure  1  -  A  schematic  of  the  geometry  considered 
in  the  derivation  of  the  reflection 
coefficient,  S?i  »  for  a  flaw.  Here, 
ft  is  the  inward  normal  to  the  flaw 
surface,  Sr  . 


For  the  situation  of  a  Rayleigh  wave  normally 
incident  on  a  half-penny  shaped  surface  crack  of 
radius  a  located  in  the  x-y  plane  (Figs.  2  and 
3),  the  reflection  coefficient  may  be  written  as 
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where  Au?  is  the  discontinuity  in  the  Rayleigh 
wave  displacement  field  across  the  crack  and  azz 
is  the  applied  stress.  The  integral  is  now  taken 
over  just  the  semi-circular  area,$  .  The  theory  we 
have  developed  for  this  configuration  is  strictly 
valid  only  in  the  low  frequency  regime  (i.e.  we 
require  the  maximum  depth  to  which  the  crack  extends 
below  the  sample  surface  to  be  much  less  than  an 
acoustic  wavelength).  To  take  into  account  the 
effect  of  imaging  at  the  surface  in  increasing  the 
value  of  the  stress  intensity  factor  near  the  sur¬ 
face,  as  well  as  the  variation  with  depth  of  the 
Rayleigh  wave  stress  fields,  we  use  the  results  of 
Budiansky  and  O'Connell7  to  write  the  surface 
integral  in  Eq.  (2)  in  terms  of  a  contour  integral 
of  the  square  of  the  mode  I  stress  intensity  factor, 
K\  ,  around  the  crack  circumference  C  (Fig.  3) 
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Here,  v  is  Poisson's  ratio  and  E  is  Young's 
modulus.  The  angular  dependence  of  K\  caused  by 
the  surface  imaging  forces  may  be  approximated  from 
the  results  of  Smith,  Emery,  and  Kobayashi.^  Smith 
et  al.  considered  the  case  of  a  half-penny  shaped 
surface  crack  in  a  beam  of  thickness  2c  subject 
to  a  bending  load.  The  applied  stress  then  takes 
the  linear  form 


aZ2(y)  =  A(1  -  y/c) 


(4) 


where  A  is  a  constant  and  y  is  the  distance 
from  the  sample  surface  (Fig.  2),  The  corresponding 
stress  intensity  factor  is  given  by 


Kj(e)  =  2/a7F  Aty0(e)  -  (a/c)^1(e )) 


(5) 


where  the  functions  ^q(9)  and  ^i(e)  were  numer¬ 
ically  evaluated  by  Smith  et  al.  (Fig.  4).  To  make 
use  of  these  results  for  our  case  where  the  stress 
is  due  to  the  Rayleigh  wave  and  not  a  bending  load, 
we  made  a  linear  approximation  to  the  Rayleigh  wave 
stress  field  and  so  evaluated  effective  values  for 
the  constants  A  and  c  appearing  in  Eq.  (4). 

Our  result  for  the  normalized  reflection  coefficient 
wk|Sn|  ,  for  a  semi-circular  crack,  as  a  function 
of  trie  normalized  crack  depth,  2va/\  ,  is  given  by 
the  dashed  curve  in  Fig.  5.  Here  k  =  2tt/X  is  the 
propagation  constant  of  the  Rayleigh  wave,  and  w 
is  the  width  of  the  acoustic  beam  at  the  crack. 


Figure  2  -  Scattering  geometry  for  a  Rayleigh  wave 
normally  incident  on  a  half-penny  shaped 
crack. 


HALF  -  PENNY  SHAPED  CRACK 

Figure  3  -  Schematic  of  half-penny  shaped  crack 
geometry. 


Figure  4  -  Angular  variation  of  the  functions  R 
\Jjq(6 )  and  ^(0)  (after  Smith  et  al.  ). 

EXPERIMENTAL  RESULTS 

In  our  experimental  studies,  we  used  a  commer¬ 
cial  hot-pressed  silicon  nitride  (NC- 132)  ceramic. 
The  samples  were  in  the  form  of  7.6  cm  x  2.6  cm  x 
.64  cm  plates  with  polished  surfaces  into  which 
cracks  had  been  introduced  via  a  Knoop  indentor. 


Measurements  were  made  at  a  frequency  of  about 
8-5  MHz*  ,  corresponding  to  an  acoustic  wavelength 
of  about  680  pm  .  Two  wedge  transducers  in  the 
configuration  shown  in  Fig.  6  were  used  to  excite 
and  receive  the  acoustic  signals.  The  experimental 
values  for  the  reflection  coefficient  thus  obtained 
were  then  corrected  for  diffraction  loss,  the  sur¬ 
face  wave  conversion  efficiencies  of  the  trans¬ 
ducers,  and  for  the  effect  of  the  transmitted  and 
reflected  acoustic  beams  each  being  at  an  angle  a 
from  the  crack  normal. 


Figure  5  -  Comparison  of  the  results  for  the  nor¬ 
malized  reflection  coefficient,  wk|Snj, 
as  a  function  of  the  normalized  crack 
depth  for  the  half-penny  shaped  crack 
theory  and  the  elliptical  crack  theory. 
The  normalized  crack  depth  is  defined 
for  each  theory  by  the  quantities  2Tra/X 
and  2tt( b  +  h ) /X  ,  respecti vely. 


Figure  6  -  Experimental  scattering  configuration 
used.  In  our  study,  a  *  12.5°  . 

Measurements  were  carried  out  on  both  annealed 
and  as-indented  samples.  One  sample  was  annealed 
in  air  and  another  in  a  vacuum  with  an  annealing 
period  of  about  six  hours  at  a  temperature  of 
I200°C  .  An  additional  three  samples  were  tested 
in  their  as-indented  state.  For  all  five  samples, 
we  monitored  the  change  in  the  acoustic  reflection 
coefficient  while  the  samples  were  slowly  stressed 
to  fracture  in  a  3-point  bending  jig. 

Annealed  Samples  -  We  found  that  the  cracks  in  the 
annealed  samples  are  fairly  well  modelled  by  the 
open  half-penny  shaped  crack  theory  described  in 
the  previous  section,  subject  to  the  additional 
observation  that  in  the  unstressed  state,  partial 
crack  closure  occurs  (Fig.  7).  As  shown  in  Figs. 

8  and  9,  loading  the  sample  from  zero  stress  to 
fracture  produces  a  steady  increase  in  the  acoustic 
prediction  for  the  crack  depth.  We  additionally 
observe  that  if,  after  a  period  of  loading,  the 
load  is  reversed  until  the  stress  in  the  vicinity  of 
the  crack  approaches  zero,  the  crack  size  estimates 
tend  to  return  to  very  near  the  intial  value 


observed  for  the  unstressed  sample  (Fig.  9).  Appli¬ 
cation  of  a  load  to  the  cracked  sample  thus  has  the 
effect  of  causing  the  crack  to  open,  as  indicated 
by  the  increase  in  the  acoustic  size  estimates  with 
increasing  load;  however,  subsequent  removal  of  the 
load  results  in  a  return  of  the  crack  to  its  origi¬ 
nal  partially  closed  state.  Verification  of  this 
interpretation  is  obtained  upon  inspection  of  the 
fracture  surface  (Fig.  10).  Here,  only  the  initial 
flaw  produced  by  the  Knoop  indentor  is  in  evidence, 
as  indicated  by  the  single  ring  near  the  sample 
surface.  Lastly,  we  compare  the  acoustic  estimate 
for  the  crack  depth  at  fracture  with  the  optically 
measured  crack  depth  (Table  I).  For  both  samples, 
agreement  between  these  two  values  is  quite  good, 
with  the  deviation  of  the  acoustically  predicted 
values  from  the  actual  values  being  less  than  17%  . 
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Figure  7  -  Fracture  model  for  crack  in  an  annealed 
sample.  Cg  is  the  depth  of  the  initial 
flaw  produced  by  the  Knoop  indentor.  In 
the  unstressed  state,  the  crack  is  par¬ 
tially  closed,  as  indicated. 
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Figure  8  -  Stress  versus  acoustic  estimates  of  the 
crack  depth  for  the  sample  annealed  in 
air  (sample  1).  Crack  depth  estimates 
are  obtained  using  the  half-penny  shaped 
crack  theory  (crack  depth,  C  ,  then 
equals  the  half-penny  shaped  crack 
radius,  a  ).  The  kink  observed  in  the 
curve  is  believed  to  be  an  effect  of  the 
severe  surface  oxidation  which  occurred 
during  the  annealing  process. 

As- Indented  Samples  -  The  fracture  model  for  the  cracks 
in  the  as-indented  samples  differs  considerably  from  the 
model  for  the  cracks  in  the  annealed  samples.  Firstly, 
there  is  crack  closure  at  the  samnle  surface  in  the  vic¬ 
inity  of  the  Knoop  indentation  (Fig.  11).  This  arises 
because  the  Knoop  indentation  technique  used  to 
produce  the  crack  creates  a  plastic  deformation 


region  in  the  neighborhood  of  the  indentation  in 
which  the  residual  stress  fields  have  not  been 
relieved  as  in  the  annealed  samples.  The  second 
difference  is  that  as  the  cracked  sample  is  stressed 
to  fracture,  the  crack  grows  from  an  initial  depth 
of  Cq  to  a  final  fracture  depth  of  Cm  . 
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Figure  9  -  Stress  versus  acoustic  estimates  of  the 
crack  depth  for  the  vacuum-annealed 
sample  (sample  2).  Two  unload-reload 
cycles  were  performed  on  this  sample. 
Crack  depth  estimates  were  obtained  using 
the  half-penny  shaped  crack  theory 
(crack  depth,  C  ,  equals  half-penny 
shaped  crack  radius,  a  ). 
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Figure  10  -  Typical  fracture  surface  for  crack  in 
an  annealed  sample. 

The  first  problem  that  we  encounter  is  that 
our  earlier  scattering  theory  based  on  the  model 
of  an  open  half-penny  shaped  crack  no  longer  is 
applicable  here.  We  have  not  yet  developed  a  theory 
which  accurately  models  the  surface  closure  charac¬ 
teristic  of  these  as-indented  cracks;  however,  we 
have  worked  out  a  theory  based  on  a  very  crude 
approximation  for  the  actual  crack  configuration 
which  gives  surprisingly  good  crack  size  predic¬ 
tions.  This  theory  approximates  the  actual  crack 


by  an  elliptical  crack  with  semi-major  axis  a  and 
semi-minor  axis  b  where  the  surface  closure  effect 
is  modelled  by  removing  the  elliptical  crack  center 
to  a  position  a  distance  h  below  the  sample  sur¬ 
face  (b  <  h  ;  Fig.  12).  We  can  evaluate  the 
reflection  coefficient  for  a  Rayleigh  wave  normally 
incident  on  this  type  of  crack  configuration  much 
as  we  did  earlier.  Again,  the  reflec¬ 

tion  coefficient  is  given  by  Eq.  (2).  Here,  how¬ 
ever,  we  must  use  a  generalized  form  of  Eq.  (3)  to 
relate  the  surface  integral  in  Eq.  (2)  to  the  con¬ 
tour  integral  involving  the  stress  intensity  factor, 
Kj  ,  namely 


Auz°zzdS  * 


2(1  -  v 


3E 


f p(r)Kj(0)dS.  (6) 


In  Eq.  (6),  the  quantity  p(r)  is  defined  for  a 
given  point  r  on  the  crack  circumference  C  as 
the  perpendicular  to  the  tangent  line  at  r  ,  which 
passes  through  the  origin  of  the  crack  coordinate 
system.  Shah  and  Kobayashi9  have  evaluated  the 
angular  dependence  of  the  stress  intensity  factor 
for  this  type  of  crack  subject  to  a  bending  load. 
The  form  for  the  applied  stress  is  again  linear 


a"z(y)  = 


C1  +  C2 


(7) 


where  Cy  and  C2  are  constants.  In  terms  of  C\ 
and  C2  *  and  the  functions  My(0)  and  M^(9) 
numerically  evaluated  by  Shah  and  Kobayashi,  the 
stress  intensity  factor  may  be  written  in  the  form 
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v^Trb/ a  00  00 

-  (a^cos^e  +  b  sin^e)  <0,^(6) 
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+  c2ML(e) 
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k  E(k)  cos  0 


(1  +  k2)E(k)  -  (k')2K(k) 


(8) 


Here,  K(k)  and  E(k)  are  complete  elliptic 
integrals  of  the  first  and  second  kinds,  respec¬ 
tively,  and  the  quantities  k  and  k"  are  defined 
by 


k2  =  1  -  (b/a)2 


(k')2  = 


(b/a)' 


(9) 


Mj(0)  and  Mj_(e)  represent  the  stress  intensity 
magnification  factors  for  the  elliptical  crack  con¬ 
figuration  shown  in  Fig.  12  for  the  respective  cases 
of  a  uniform  applied  stress  and  a  linearly  varying 
applied  stress.  These  functions  have  been  calcu¬ 
lated  by  Shah  and  Kobayashi  for  various  values  of 
the  aspect  ratio,  b/a  ,  and  the  depth  ratio,  b/h  . 

As  in  the  actual  crack,  the  crack  width  is  roughly 
twice  the  crack  depth,  and  the  region  of  surface 
closure  is  shallow  compared  to  the  maximum  crack  depth. 
We  have  chosen  to  use  the  Shah  and  Kobayashi  results 
for  My (6 )  and  M|_(e)  given  for  the  parameter 
values  b/a  =  0.6  and  b/h  =  0.9  (Figs.  13  and 
14).  The  result  for  the  normalized  reflection 
coefficient,  wk|Sn|  ,  as  a  function  of  the  nor¬ 
malized  crack  depth  2tt ( b  +  h)/x  ,  is  then  given  by 
the  solid  curve  in  Fig.  5. 
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Figure  11  -  Fracture  model  for  an  as-indented  crack. 
Cq  is  the  depth  of  the  initial  flaw 
produced  by  the  Knoop  indentor;  Cm  is 
the  depth  to  which  the  crack  grows 
before  fracturing.  Crack  closure  occurs 
at  the  surface  in  the  vicinity  of  the 
Knoop  indentation. 


Figure  12  -  Schematic  of  elliptical  crack  theory 
geometry. 


Figure  13  -  Stress  intensity  magnification  factor 
as  a  function  of  the  angle  9  for  the 
elliptical  crack  configuration  of 
Fig.  12,  subject  to  a  uniform  applied 
stress  (after  Shah  and  Kobayashi9). 
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Figure  14  -  Stress  intensity  magnification  factor 
as  a  function  of  the  angle  6  for  the 
elliptical  crack  configuration  of 
Fig.  12,  subject  to  a  linearly  varying 
applied  stress  (after  Shah  and 
Kobayashi9) . 

We  used  this  crude  theoretical  model  to  give 
crack  size  estimates  for  the  as-indented  samples. 

As  with  the  cracks  in  the  annealed  samples,  there 
is  a  steady  increase  in  the  crack  depth  estimates 
with  increasing  load  (Figs.  15,  16,  and  17).  How¬ 
ever,  with  these  cracks,  crack  growth  occurs,  as 
indicated  by  the  failure  of  the  cracks  to  return  to 
their  original  sizes  upon  subsequent  unloading.  In 
addition,  a  small  decrease  in  the  crack  size  esti¬ 
mates  is  observed  when  the  load  is  relieved,  thus 
indicating  a  slight  tendency  of  the  cracks  towards 
partial  closure  at  the  lower  boundary  upon  unload¬ 
ing.  Continuous  growth  of  the  crack  until  fracture 
is  evidenced  in  Fig.  16  where  two  unload-reload 
cycles  were  performed.  In  each  cycle,  it  is  appar¬ 
ent  that  the  crack  never  fully  returns  to  the  size 
it  had  previously  when  the  stress  is  reduced  to  near 
zero  values. 
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Figure  15  -  Stress  versus  acoustic  estimates  of  the 
crack  depth  for  as-indented  sample  4. 
One  unload-reload  cycle  was  performed. 
The  crack  depth,  C  ,  equals  the 
quantity  b  +  h  (see  Fig.  12). 
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Figure  16  -  Stress  versus  acoustic  estimates  of 
the  crack  depth  for  as-indented 
sample  5.  Two  unload-reload  cycles 
were  performed.  The  crack  depth,  C  , 
equals  the  quantity  b  +  h  (see 
Fig.  12). 
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Figure  17  -  Stress  versus  acoustic  estimates  of  the 
crack  depth  for  as-indented  sample  6. 

One  unload-reload  cycle,  where  the 
sample  was  fully  unloaded,  was  per¬ 
formed.  The  crack  depth,  C  ,  equals 
the  quantity  b  +  h  (see  Fig.  12). 

Examination  of  a  typical  fracture  surface  for 
these  samples  reveals  the  presence  of  two  rings 
near  the  sample  surface  (Fig.  18).  The  inner  ring 
represents  the  extent  of  the  initial  flaw  produced 
at  the  time  of  indentation.  The  outer  ring  indi¬ 
cates  the  extent  to  which  the  crack  grew  before 
fracture  occurred.  A  comparison  of  the  acoustic 
estimates  for  the  initial  crack  depth,  Cq  ,  and 
the  depth  at  fracture,  Cm  ,  with  the  corresponding 
values  measured  optically  is  given  in  Table  II. 
Considering  the  crudeness  of  the  model  used  to 
obtain  the  acoustic  estimates,  the  agreement 
between  these  values  is  remarkably  good. 

Lastly,  we  compare  the  fracture  stresses  of 
the  five  samples.  From  Tables  I  and  II,  we  note 
that  the  crack  depths  at  fracture  for  the  annealed 
samples  tend  to  be  much  less  than  those  for  the 
as-indented  samples.  Thus,  we  expect  the  annealed 
samples  to  fracture  at  higher  stress  values  than 
the  as-indented  samples.  In  Table  III,  we  show  that 
this  is  indeed  the  case. 


Figure  18  -  Typical  fracture  surface  for  crack  in 
an  as-indented  sample. 

CONCLUSION 

We  conclude  that  the  acoustic  technique  is  a 
powerful  method  for  observing  fracture  of  surface 
cracks  in  ceramics.  Major  differences  in  the 
fracture  behavior  of  cracks  in  annealed  and 
unannealed  samples  have  been  documented. 

It  is  apparent  that  there  is  a  need  to  improve 
the  acoustic  scattering  formulae  and  the  estimates 
of  the  stress  intensity  factor  for  cracks  closed  at 
the  sample  surface.  A  theory  by  which  to  accurately 
predict  how  a  partially  closed  annealed  crack  opens 
under  applied  stress  would  also  be  of  considerable 
interest.  We  are  currently  developing  a  variational 
method  to  deal  with  the  problems  encountered  in  the 
scattering  theory  and  are  considering  simple 
theories  to  deal  with  the  latter  problem. 
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Acoustic 

Actual 

Sample 

Cm  <ym) 

Cm  (ljm> 

Deviation 

1 

134 

115 

17" 

2 

*25 

115 

7% 

7.  8.  Budiansky  and  R.J.  O'Connell,  "Elastic 
Moduli  of  a  Cracked  Solid,"  Int.  J.  Solid 
Structures,  vol.  12,  pp.  81-97,  Pergamon  Press, 
1976. 

8.  F.W.  Smith,  A.F.  Emery,  and  A.S.  Kobayashi, 
"Stress  Intensity  Factors  for  Semicircular 
Cracks:  Part  2,  Semi -Infinite  Solid,"  J.  Appl . 
Mech.,  vol.  34,  no.  4,  pp.  953-959,  December 
1967. 


9.  R.C.  Shah  and  A.S.  Kobayashi,  "On  the  Surface 
Flaw  Problem,"  The  Surface  Crack  Physical 
Problems  and  Computational  Solutions, 

J.L.  Swedlow,  ed.  pp.  79-124;  presented  at  the 
Winter  Annual  Meeting  of  the  ASME,  New  York, 

New  York,  November  1972. 

10.  G.S.  Kino,  "Variational  Formulae  for  Scattering 
of  Acoustic  Waves  by  Flaws  and  for  Acoustic 
Wave  Propagation,"  Ginzton  Laboratory  Report 
No.  2634,  Stanford  University,  November  1976. 

11.  G.S.  Kino,  "Variational  Methods  for  Calculating 
J  and  M  Integrals  and  Acoustic  Scattering 
From  Cracks,"  Materials  Research  Council  Report, 
July  1980. 


Sample  Fracture  Stress  (MPa) 


Annealed:  1  332 

2  340 

As-Indented:  4  238 

5  239 

6  242 


References 

l;  J.  Tien,  B.  Khuri-Yakub,  and  G.S.  Kino,  "Acous¬ 
tic  Surface  Wave  Probing  of  Ceramics,"  Proc. 
ARPA/AFML  Rev.  of  Progress  in  Quantitative  NDE, 
La  Jolla,  California,  July  1979. 

2.  O.B.  Marshall  and  B.R.  Lawn,  "Residual  Stress 
Effects  in  Sharp  Contact  Cracking:  Part  1, 
Indentation  Fracture  Mechanics,"  J.  Matls. 
Science,  vol.  14,  pp.  2001-2012,  1979. 

3.  D.B.  Marshall,  B.R.  Lawn,  and  P.  Chantikul, 
"Residual  Stress  Effects  in  Sharp  Contact  Crack¬ 
ing:  Part  2,  Strength  Degradation,"  J.  Matls. 
Science,  vol.  14,  pp.  2225-2235,  1979. 

4.  B.R.  Lawn  and  D.B.  Marshall,  "Residual  Stress 
Effects  in  Failure  From  Flaws,"  J.  Am.  Ceramic 
Soc.,  vol.  62,  no.  1-2,  pp.  106-108,  Jan-Feb 
1979. 

5.  G.S.  Kino,  "The  Application  of  Reciprocity 
Theory  to  Scattering  of  Acoustic  Waves  by 
Flaws,"  J.  Appl.  Phys.,  vol.  49,  no.  6, 
pp.  3190-3199,  June  1978. 

6.  B.A.  Auld,  "General  Electromechanical  Recipro¬ 
city  Relations  Applied  to  tne  Calculation  of 
Elastic  Wave  Scattering  Coefficients,"  Wave 
Motion,  vol.  1,  no.  1,  pp.  3-10,  January  1979. 


cfi> 


Contract  N00014-78-C-0283 
Contract  Period:  1  March  1980  -  28  February  1981 
Amount  of  Contract:  $200,812 


Scientific  Officer: 

Dr.  Arthur  M.  Diness 
Director,  Metallurgy  and 
Ceramics  Program 
Materials  Sciences  Division 
Office  of  Naval  Research 
Department  of  the  Navy 
800  North  Quincy  Street 
Arlington,  Virginia  22217 


G.  L.  Report  No.  3256 
April  1981 

Edward  L.  Ginzton  Labroratory 
W.  W.  Hansen  Laboratories  of  Physics 
Stanford  University 
Stanford,  California 


6  12  076 


>o 


MEASUREMENTS  OF  SURFACE  DEFECTS  IN  CERAMICS 


Introduction 

During  the  past  year  we  have  continued  our  acoustic  surface  wave 
measurements  of  the  scattering  coefficients  of  surface  cracks  in  ceramics. 

The  aim  is  to  measure  the  size  of  the  surface  cracks  and  predict  the  failure 
stress  of  parts  containing  them. 

The  research  has  been  extremely  successful ;  we  can  now  make  accurate 
predictions  of  failure  stress  of  unannealed  samples  containing  penney-shaped 
cracks  in  the  100-300  um  range.  By  changing  the  frequency,  the  basic 
technique  should,  in  fact,  be  usable  over  a  far  greater  range  in  both 
directions(directions?)  of  crack  sizes.  Not  only  can  we  predict  the  failure 
stress  when  the  applied  stress  is  perpendicular  to  the  direction  of  the  crack, 
but  we  can  also  make  predictions  accurate  to  within  20%  of  the  failure 
stress  of  cracks  at  angles  up  to  45%  to  the  applied  stress,  provided  we  make 
the  acoustic  wave  measurement  in  the  direction  of  the  applied  stress. 

We  have  studied  the  effects  of  residual  stress  on  crack  size  estimation, 
fracture  stress  and  slow  crack  growth,  and  now  have  a  very  good  understanding 
of  the  effects  of  plasticity  and  residual  stress  in  ceramics.  The  acoustic 
measurements  have  enabled  us  to  watch  a  crack  grow  and  to  examine  the  effects 
of  crack  closure.  Thus,  they  have  provided  an  important  new  experimental  tool 
for  fundamental  studies  of  crack  growth  in  ceramics.  We  are  presently 
extending  our  study  to  evaluate  the  effects  of  surface  roughness  and  machine 
damage  on  crack  size  estimation  and  fracture. 


Research  Progress 

The  acoustic  surface  wave  measurement  of  the  surface  cracks  is  carried 
out  at  a  frequency  of  8  MHz  ,  and  the  cracks  range  in  size  from  50  to 
300  um  .  Thus,  we  operate  in  the  region  of  ka  =  .5  to  ka  -  3  .  We  have 
already  reported  on  the  theory  developed  to  measure  the  size  of  surface  cracks 
from  a  scattering  measurement  in  the  low  ka  regime.  However,  our  initial 
theory  did  not  take  into  account  the  effect  of  residual  stress  on  crack  sizing 
and  fracture  prediction.  This  deficiency  in  the  theory  was  evidenced 
experimentally  by  undersizing  the  surface  cracks  by  a  factor  of  2  or  more. 

At  the  end  of  the  last  reporting  period,  we  realized  that  due  to  the  residual 
stress  (introduced  by  the  Knoop  indentor),  the  crack  tip  is  closed  at  the 
surface,  and  the  reflection  coefficient  of  the  crack  is  greatly  reduced.  We 
have  now  developed  a  model  for  such  cracks,  and  can  theoretically  predict  the 
behavior  of  surface  cracks  in  the  presence  or  absence  of  residual  stress.  We 
have  confirmed  this  theory  experimentally,  and  can  now  size  cracks  and  predict 
the  failure  of  parts  containing  cracks  with  an  accuracy  of  less  than  20%  . 

Appendix  A  is  a  paper  that  will  be  published  in  the  ARPA/AFML  Review  of 
Progress  in  Quantitative  NDE,  1980,  in  which  we  describe  the  theory  and 
experimental  results  on  the  measurement  of  surface  cracks  in  silicon 
nitride.  In  summary,  the  Knoop  indentor  introduces  residual  stress  that 
closes  the  crack  at  its  surface.  Thus,  the  crack  displacement  is  reduced,  and 
the  reflection  coefficient  from  the  crack  is  reduced  as  compared  to  a  crack 
with  open  crack  tips.  As  a  first  order  approximation,  we  consider  the  crack 
to  be  elliptical  in  shape  lying  just  below  the  surface  (10%  from  total  crack 
depth).  We  calculate  the  reflection  coefficient  of  an  elliptical  embedded 


crack,  and  we  use  this  to  predict  crack  size.  Measurements  of  crack  size  were 
made  on  several  samples  with  and  without  (annealed)  residual  stress.  We 
measured  the  reflection  coefficients  of  the  surface  cracks  as  they  were 
stressed  to  fracture  in  a  three-point  bending  jig. 

Our  experiments  indicate  that  annealed  samples  behave  like  surface 
breaking  surface  cracks,  with  open-tips  (no  residual  stress),  and  that  no  slow 
crack  growth  is  present,  as  would  be  expected  in  brittle  materials.  The  crack 
sizes  and  fracture  stresses  were  predicted  with  an  accuracy  of  better  than 
20%  (Table  I,  Appendix  A). 

In  unannealed  samples,  we  found  evidence  of  slow  crack  growth.  Slow 
crack  growth  is  observed  as  shown  in  Figs.  15,  16,  and  17  of  Appendix  A,  where 
the  crack  does  not  go  back  to  its  initial  size  when  unloading  the  bending 
stress  in  the  three-point  bending  jig.  These  cracks  behaved  like  buried 
elliptical  cracks  just  below  the  surface.  The  estimate  of  crack  size  and 
fracture  stress  was  again  predicted  with  an  accuracy  of  better  than  20%  .  We 
have  thus  a  first  order  model  that  predicts  the  behavior  and  effect  of  surface 
cracks  when  their  annealing  history  is  known. 

A  series  of  experiments  were  carried  out  to  measure  the  reflection 
coefficients  of  surface  cracks  that  were  oriented  at  an  angle  to  the  long  axis 
of  the  rectangular  test  bars.  The  angle  between  the  long  axis  of  the  test 
bars  and  the  cracks  varied  from  90  to  45u  .  The  samples  were  tested  with  a 
fixed  angle  of  incidence  of  the  acoustic  waves  to  the  long  axis  of  the  bars, 
and,  the  samples  were  broken  in  the  same  three-point  bending  jig.  No  attempt 
was  made  to  predict  the  crack  size  or  its  growth  as  the  samples  were  stressed 
to  fracture.  Only  the  reflection  coefficient  was  measured.  The  crack  sizes 


were  not  estimated  because  the  cracks  will  grow  in  a  direction  perpendicular 
to  the  direction  of  the  applied  stress,  thus  giving  skew-shaped,  not  planar 
cracks.  The  reflection  coefficient  from  the  cracks  was  found  on  the 
assumption  that  an  equivalent  crack,  normal  to  the  direction  of  incidence  of 
the  acoustic  wave,  existed,  and  the  failure  stress  was  predicted  on  the  basis 
of  the  same  assumption.  The  actual  fracture  stress  was  predicted  with  an 
accuracy  of  better  than  10%  for  the  10  cracks  tested.  This  result  is 
extremely  encouraging.  We  had  indeed  expected  it  because  the  long  wavelength 
scattering  theory  we  used  indicated  that  the  fields  of  the  surface  waves  are 
perturbed  in  the  same  manner  as  the  field  of  the  static  stress  applied  to 
fracture. 

At  the  present  time,  we  are  carrying  out  a  series  of  experiments  designed 
to  evaluate  the  effect  of  machining  and  surface  roughness  on  the  strength  of 
silicon  nitride  in  the  presence  or  absence  of  surface  cracks,  and  in  the 
presence  or  absence  of  residual  stress.  The  samples  are  being  machined  and 
should  be  tested  shortly. 

Part  of  the  error  in  estimating  crack  sizes  from  the  reflection 
coefficient  measurement  is  due  to  the  model  used  for  the  crack  with  closed 
tips  (buried  ellipse).  We  have  started  working  on  a  variational  method  to 
predict  the  reflection  coefficient  of  cracks  with  odd  shapes  and  mixed 
boundary  conditions  (closed  or  open).  Appendix  B  is  a  pre-print  summarizing 
the  variational  technique.  So  far  we  have  tried  the  theory  on  an  infinitely 
long  strip  crack  with  a  known  solution  to  the  crack  displacement.  By  using  a 
trial  guess  at  the  displacement,  we  predict  a  reflection  coefficient  that  is 
10%  lower  than  the  exact  solution,  which  results  in  an  error  of  less  than 


-4- 


1 


2%  in  the  value  of  the  stress  intensity  factor.  The  variational  theory  is 
being  programmed  so  that  an  optimum  computer  solution  can  be  generated.  We 
also  plan  to  expand  the  theory  to  evaluate  the  reflection  coefficient  of 
two-dimensional  cracks  of  odd  shapes. 
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APPENDIX  B 


VARIATIONAL  METHODS  FOR  CALCULATING  J  &  M 
INTEGRALS  AND  ACOUSTIC  SCATTERING  FROM  CRACKS 
G.S.  KINO 


ABSTRACT 


It  is  difficult  to  calculate  acoustic  scattering 
coefficients,  stress  intensity  factor  and  energy  integrals  of 
cracks  of  arbritrary  shape  and  with  non  uniform  applied  stress 
fields.  We  derive  here  a  variational  technique  for  this  pur¬ 
pose,  which  uses  a  trial  function  for  the  crack  opening  dis¬ 
placement  AU.  We  show  that  the  method  is  easily  adaptable  to 
the  use  of  fourier  transform  techniques  and  hence  can  make  use 
of  the  FFT  (fast  fourier  transform)  numerical  methods  commonly 
employed  in  other  fields.  We  give  a  very  simple  analytic  ex¬ 
ample  of  a  plane  crack  in  mode  I  plane  strain.  Using  a  crude 
triangular  trial  function  for  the  crack  opening  displacement 
we  obtain  interaction  energy  (M  integral)  15%  lower  than  the 
true  value,  with  a  corresponding  stress  intensity  factor  8% 
lower  than  the  true  value.  We  suggest  the  use  of  Rayleigh-Ritz 
techniques  to  obtain  more  accurate  results. 
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VARIATIONAL  METHODS  FOR  CALCULATING  J  &  M 
INTEGRALS  AND  ACOUSTIC  SCATTERING  FROM  CRACKS 

G.S.  KINO 

INTRODUCTION 

It  is  normally  difficult  to  calculate  the  stress  intensity 
factor  of  flat  crack  of  arbitrary  shape.  Similarly,  the  J  &  M 
integrals  or  the  closely  related  scattering  coefficient  of  an 
acoustic  wave  reflected  from  a  crack  can  only  be  calculated 
easily  for  certain  simple  crack  shapes  (ellipses  or  strips) .  A 
further  difficult  problem  is  to  determine  these  quantities  when 
the  applied  stress  is  non-uniform  over  the  area  of  the  crack. 

At  the  1975  MRC  conference,  the  author  derived  an  acoustic 
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wave  scattering  theory  '  '  which  could  be  stated  in  variational 

form.  In  essence  the  use  of  such  variational  theories  could  make 
it  possible  to  employ  a  trial  function  AU  for  the  crack  opening 
displacement.  A  more  accurate  guess  is  then  obtained  for  the 
energy  related  quantity  of  interest  such  as  the  M  integral  or 
acoustic  wave  scattering  coefficient.  The  variational  method 
has  the  advantage  of  making  it  possible  to  use  physical  intuition 
in  the  choice  of  the  trial  function  and  Rayleigh-Ritz  techniques 
to  improve  the  trial  function.  Furthermore,  it  lends  itself, 
as  we  shall  see,  to  the  use  of  FFT  (Fast  Fourier  Transform  Tech¬ 
niques)  and  makes  it  possible  to  employ  many  of  the  powerful  methods 
that  have  been  used  earlier  in  electromagnetic  scattering  theories. 


Such  techniques  do  not  involve  the  use  of  singular  integral 
equations  nor  of  Wiener  Hopf  methods.  So,  mathematically  they 
are  relatively  simple  and  suitable  for  numerical  computation,  and 

are  not  restricted  only  to  elastostatic  solutions.  The  theory 
may  therefore  be  useful,  ultimately,  in  applications  to  acoustic 
wave  scattering  theory  at  arbitrary  frequencies. 

In  this  paper  we  have  re-derived  the  variational 
theory  by  a  simpler,  hopefully,  more  easily  understood  method. 

We  show  that  it  may  be  generalized  to  determine  not  only 
the  scattering  coefficient  of  acoustic  waves,  but  also  the  in¬ 
teraction  energy  of  a  crack  and  its  M  integral. 

We  then  go  on  to  demonstrate  the  use  of  fourier  transform 
techniques  on  a  simple  problem,  that  of  a  plane  crack  in  static 
Mode  I  plane  strain.  We  show  that  the  use  of  a  relatively  crude 
triangular  trial  function  for  the  crack  opening  displacement  AU  = 
a- | x | ,  instead  of  the  true  function  AU  =  (a  -x  )  where  2a  is 
the  crack  length  and  x  the  distance  along  its  surface,  gives  an 
interaction  energy  and  stress  intensity  factor,  15%  and  8%  lower 
than  the  true  values  respectively. 


SCATTERING  COEFFICIENT  AND  ENERGY  INTEGRALS 


We  employ  the  Rayleigh-Betti  reciprocity  theorem  to  determine 

the  energy  change  due  to  the  presence  of  a  flaw.  We  consider  the 

material  to  be  stressed  in  a  testing  machine.  Suppose  that  without 

a  flaw  present  the  displacement,  stress  and  strain  at  any  point 

a  A  A 

in  the  material  are  Uv,  a.  .,  and  e. •,  respectively.  With  the 

i  i3  ij  1 

flaw  present  the  displacement  and  stress  at  IL  ,  0  ^y  e±j'  respectively . 

By  the  usef  of  the  Rayleigh-Betti  reciprocity  theorem  in  the 
material  region  outside  the  flaw,  it  follows  that: 


(0tjuj  -  nids  “ 
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0  .  .U  . 
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n .  ds 
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(1) 


Where  S  is  the  outer  surface  of  material,  S?  the  surface  of  the 
flaw,  and  n^  the  outward  normal  from  these  surfaces. 

We  note  that  the  only  contribution  to  the  integral  S  comes 
from  the  two  ends  of  the  material  where,  we  assume  it  is  supported 
in  rigid  clamps.  (o^.n^  =  0  at  a  free  surface) .  Then  at  the 

A 

clamps  =  U ^ .  So  it  will  be  seen  that  the  interaction  energy 
with  a  flaw  is  W  where: 
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When  the  flaw  is  a  crack  o^jn^  =  0  on  the  crack  surface.  So  it 
follows  that  the  crack  energy  is 
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Where  AUj  is  the  crack  opening  displacement,  is  the  outward  normal 
on  one  side  of  the  crack.  The  energy  of  a  more  general  flaw  is 
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These  results  agree  with  those  of  Eshelby 

It  follows  from  Budiansky  and  O'Connell5  that  the  interaction 
energy  of  a  crack  is 
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where  p  and  l  are  based  on  the  crack  based  coordinates  shown  in 
Fig  (1) .  E  is  Young's  modulus,  v  is  Poisson's  ratio  and  K ^ ,  K_^, 
and  are  the  stress  intensity  factors  for  plane  strain  and 

antiplane  shear.5  In  turn,  Rice's  J  integral6  of  fracture  mechanics 
is  defined  as 


j  =  1=21 
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while  the  M  integral  of  Knowles  and  Sternberg  is  defined  by  the 
relation 


M  = 


Jd  l 


(7) 


It  therefore  follows  that 

W  =  M/3  (8) 

We  conclude  that  if  we  can  evaluate  w  we  can,  in  turn  evaluate 
M.  Furthermore,  it  was  shown  by  Kino  that  when  an  acoustic  wave  is 
scattered  from  a  crack,  the  back  scatter  reflection  coefficient 
of  the  wave  of  frequency  a»  is  ' 
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where,  now,  U.,  a- •  are  the  values  of  U,o. ■  at  the  crack  when 
3  13  3  13 

there  is  unit  power  emitted  by  a  transducer.  So  S^  is  just  j  to  W/4 
with  an  applied  stress  equivalent  to  that  of  unit  power  in  the 
incident  acoustic  wave. 

Therefore,  if  we  can  evaluate  W  correctly,  we  can  determine 
M  and  S^,  parameters  of  importance  in  fracture  mechanics  or  in 
the  determination  of  acoustic  wave  scattering  from  a  crack.  For 

A 

an  exact  calculation  we  need  to  know  the  applied  field  a ^  .  and 
the  crack  opening  displacement  AU^  in  the  pressure  of  the  applied 
field.  Thus  it  is  necessary  to  solve  for  AU.  when 


AU  .  =  0 
3 


outside  the  crack 


a  .  .n  .  =  0 
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at  the  crack  surface 


and 


o .  .  is  known. 
13 


VARIATIONAL  THEORY 


We  have  found  that  the  energy  stored  in  a  crack  is 


W  =  (oA .U .  -  a.  .UA  )n.ds 
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Suppose  we  use  a  trial  function  for  the  fields  within  the  interior 

in  which  U.,a. .  etc.  differ  to  1st  order  error  from  the  correct 
i  ID 

quantities  tL,a^  .  Then  if  we  consider  the  quantity,  called  the 
reaction,  by  Rumsey^ 


AW  =  (o.  .U  -  a  .U .  )  n.ds 
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taken  around  the  outer  periphery,  we  see  that  when  a.  .  =  a.  ., 
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Suppose  now  we  require  AW  to  differ  from  zero  only  to  2nd 

order  in  error.  Then  if  is  the  applied  displacement  at  the 

grips  U .  =  U°  =  UA. 

D  D  D 
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It  therefore  follows  that  if  Aw  =  0  to  2nd  order  in  error, 
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is  a  formula  for  w»  a  correct  to  2nd  order  in  error  when  trial 
functions  Uj,o^  are  employed. 

But  it  follows  from  the  Betti-Rayleigh  reciprocity  theorem 


that 
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at  the  flaw. 


If  the  flaw  is  a  crack  .n  =  0  at  the  crack  surface  and  so 
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Aw  becomes 


AW  =  -  (AU°.  a.  .  )n  .  ds  =  0 
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is  the  requirement  for  zero  1st  order  error  in  w  ,  where  n  j_  is  the 
outward  normal  from  one  surface  of  the  crack.  But  0.^  ^  is  only 
non  zero  at  the  crack  to  1st  order  in  error.  Hence  we  can  replace 
AU?  by  AUj  and  only  contribute  terms  in  error  to  2nd  order  in  Aw.  So 


AW  =  -  (AU  .0 .  . ) n.  ds 
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is  an  equally  valid  assumption. 

It  follows  that  Eq.  (16)  is  the  requirement  for  our  variational 
principle.  It  requires  continuity  of 


(ai jAUj>nids 


through  the  crack  instead  of  the  far  more  stringent  boundary 


condition  :o^.n^ 
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Mow  consider  W. 


When  we  use  a  trial  function  AU^,  it  follows 


that 
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Where  the  total  field  is 

A  ,  s 
a .  .  =a  .  .  +  a  .  . 
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So  a  is  the  stress  field  due  to  the  displacement  AU . . 
13  1 


Because 


(UijAUj) n^ds  =  0 


we  see  that 
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But  as  we  do  not  know  the  amplitude  AU^  and  is  proportional 

to  AU  j  ,  we  could  equally  well  replace  AU_.  by  AAU^  ,  where  A  is  a 
proportionality  factor.  In  this  case  it  follows  that 


W  =  -  a 
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However  we  could  equally  well  have  written 
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In  this  case 
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It  follows  that  as  Eqs .  (20)  and  (22)  must  be  identical  to  2nd 


order  hence 
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This  is  the  variational  formula  required. 

We  can  prove  that  this  formula  is  variational  by  writing 
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where  x»  x'are  coordinates  on  the  crack  surface,  and  the  Green's 
function  G^(x,x")  is  normally  symmetric  in  the  x  and  x'  cooic.  .es 
It  follows  from  Eqs. (25)  and  (26)  that 
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Substitution  from  Eqs.  (21)  and  (26)  in  Eq.  (27)  leads  to  the 
result  6w  =  0.  So  the  formula  of  Eq.  (25)  is  variation. 


The  theory  can  be  generalized  fairly  easily  using  the 
methods  of  reference  1  to  show  that  the  variational  formula 
for  scattering  from  a  crack  by  a  signal  incident  from  transducer 
1  and  received  on  transducer  2  is 
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Where  the  supercripts  (1)  and  (2)  refer  to  the  fields  associated 
with  incident  waves  from  transducers  1  and  2 / respectively . 


FOURIER  ANALYSIS  AND  VARIATIONAL  TECHNIQUES  FOR  A  CRACK 


The  variational  methods  lend  themselves  well  to  the  use 
of  Fourier  analysis  techniques,  and  in  particular  the  application 
of  the  Fast  Fourier  transform  numerical  computer  routines  which 
are  now  widely  available. 

Our  basic  approach  will  be  to  use  a  trial  function  for 

AU  on  the  crack  surface  z  *  0.  Then  on  the  plane  z  =  0  we  can 
2 

write 
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It  is  relatively  easy  to  write  down  the  Fourier  transform  of 
other  quantities  Ux,Uy/Uz,  ozz>oxx,o  .  Then  at  any  point 


(31) 


in  space  we  can  use  the  plane  wave  angular  spectrum  analysis  to 
find  the  connections  between  these  various  quantities. 

If  for  instance,  a  (x,y,o)  has  the  Fourier  transform 

zz 

it  follows  from  Rayleigh's  theorem  that 

fa  s  f  * 

oV,  AU  dxdy  =  o (ct,3)AU  (a,S)dadB 

Thus  it  is  easy  to  evaluate  the  terms  we  need  in  the 
variational  theory  by  relatively  simple  integrations,  and  to 
treat  flat  cracks  of  arbitrary  shape. 

We  shall  use  as  an  illustration  of  the  variational 
technique  derived  here  an  example  calculated  for  scattering  in 
unpublished  work  by  Kino  and  Tien,  using  the  earlier  1975  formu¬ 
lation  of  the  variational  theory.  We  take  the  crack  to  be  located 
in  the  z=0  plane,  to  be  of  length  2a  in  the  direction  and  infinitely 
long  in  the  y  direction,  as  illustrated  in  Fig.  2. 

We  shall  treat  only  the  static  case  of  a  ^laiHe , 
crack  with  mode  I  excitation.  It  will  be  clear  from  the  treatment 
how  to  deal  with  more  complicated  situations  with  time  varying 
excitations  and  more  complicated  shapes.  We  shall  compare  our 
results,  using  simple  trial  functions,  with  the  exact  solution. 

We  write  the  displacement  in  the  form 

U  =7$  +  (32) 

where 
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An  illustration  of  the  plane  crack  used  in  the  theory. 


with 
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for  waves  whose  components  vary  as  exp  jut.  The  mass  density 

of  the  material  isp  and  the  Lam4  constants  Xand  u . 

mo 

The  potentials  on  each  side  of  the  crack  are  taken  to 
be  4>+  (x,  z)  ;  ^+(x,y)  ;  and  <j>  (x,z),  \p  (x,z),  respectively. 

We  can  carry  out  Fourier  transforms  in  the  x  direction 
of  these  quantities  and  write 


+  /  . 

4>  (x,z)  = 

f 

A  (a) 

e-jax 

e'i!lz  da 

(35) 

4j+(x,z)  = 

B  (a) 

e'jax 

e-j8s2  da 

(36) 

<f>“(x,z)  = 

C  (a) 

e~jax 

e?^i,z  da 

(37) 

4 ~ (x,z)  = 

D  (a) 

e“jaX 

da 

(38) 

where 

2  =  k2  (39) 


2  2 
a  +es 


(40) 


It  will  be  seen  that 


A(a)  =  J- 


<J>+(x,o)  e^ax  dx 


(41) 


B  (a)  -  57 


^+(x,o)  e^aX  dx 


(42) 


C  (a)  =  ~ 

2  TT  J 


<J>"(x, o)  eJotX  dx 


(43) 


°(a)  =  f- 


ip~  (x,o)  e-Iax  dx 


(44) 


In  order  to  find  a j,  we  take  the  boundary  conditions  at  the  crack 
to  be 


AU  =  U  (x,o)  -  U  (x , o ) 

Z  w  z 


(45) 


A Ux  =  Ux(x,o)  -  Ux(xfo)  =  0 


(46) 


°zz{x'o)  =°zz(x'o)  =°zz 


(47) 


°xz(X'o)  =Gxz(x'o) 


(48) 


The  use  of  Eq.  (45)  with  Eqs .  (41),  (42),  and  (32)  yields  the 

result 

Jax  (49) 


8l(A  +  C)  +a(B  -  D)  =  L- 


AU  eJ~'*  dx 
z 


Similarly  from  Eq.  (46)  it  follows  that 


0 


(50) 


a  (A  -  C)  -6  (B  +  D)  = 


It  may  similarly  be  shown  from  Eq.  (48)  that 


2aQl  (A  +  C)  +  (a2  -  Bg)  (B  -  D)  =  0 


(51) 


Where  we  have  used  the  relation 


l+2u  =  _£  = 


2  2 
a  +B 
_ s_ 

2  .  n  2 

“  +62 


(52) 


Continuity  of  o22  Q Eq.  (470  yields  the  result 


(B?  ~  a2) (A  -  C)  +  2ag  (B  +  D)  =  0 
s  s 


(53) 


It  follows  from  Eqs.  (50)  and  (53)  that 


A  =  C 


(54) 


D  =  -B 


(55) 


So  from  Eqs.  (51) ,  (54)  and  (55)  we  see  that 


2aB , 


B  =  - 


2  a  2 

a  -8S 


(56) 


It  therefore  follows,  after  substitution  of  Eqs.  (54),  (55) 

and  (56)  in  Eq.  (49)  that 


**  H 


+a' > j 


with 


AU2  =  -2j 


6 1  -  a2  Ae^ac 


1  „2  2 

Bs  -a 


and 


CT  =  ~V 
zz 


p 

(g  -a  )  +  4a  g  g,,  -] 

A 

J  L. 

c  2  2  \ 
6  -a  -1 

-lax  , 
e  J  da 


It  follows  that 


a  A  U  dx 
zz  z 


=  +2ju 


|A  0‘ 


(B2+a2)  [  (g2-a2)2+4a2gsS j 


2  2  2 
(8g  -  o' V 


da 


We  note  that  g  and  g  must  always  be  positive,  i-e.  we  take 

s  t 


their  magnitude-  So  it  is  convenient  to  write 


o^AUzdx  =  4j, 


,2  <4+ 


a2)  [  (  gj-a2)  +  4a2g  gn  J  e 
s  s  a,  u  a 


i 


.  2  2,2 

(  Ss  -  a  ) 


If  A (a)  is  known,  this  integral  can  be  evaluated  numerically. 
We  put 


a  2  2 
gs  -  a 


60(B2  +a2) 


(58) 


(59) 


(60) 


(61) 


A (a )  =  R(a) 


(62 


Where 


R  =  AU  eDax  dx 

4tt  z 


So  it  follows  that 


a  AU  dx 
zz  z 


=  8jyir 


Qg^-a2)2  +  4a2S  S,  J 


SjteJ  +  a2) 


We  shall  now  evaluate  the  static  case  for  which  k  -*■  0 , 

s 

0.  In  this  case  the  integral  reduces  to  the  form 


a  AU  dx  = 
zz  z 


-  16yrr 


‘a  1- 


r"  AU  dx  = 
zz  z 


4Ett  f  ,_2, 

- 5—  I R  J  a  da 

(1-0  J 

0 


where  E  is  Young's  modulus,  v  is  Poisson's  ratio. 
Thus  it  follows  that 


[f°j?  ]  2 

flR!2a  da 

1-v  J 


Where  R  is  given  by  Eq.  (63) . 

We  now  use  the  variational  formula  of  Eq.  (67)  to  evaluate 

W.  We  assume  a  =  constant,  initially. 

zz  J 


Example  I;  Exact  Solution 
AUz  =  KU2-*2,1''2 


x  <  a 


(68) 


AU  =  0 
z 


x  >  a 


In  this  case  it  can  be  shown  that 


=  jK(a/4a)  J^aa) 


(69) 


Where  J^(x)  is  a  1st  order  Bessel  function  of  the  1st  kind. 
So  it  follows  that 


W  = 


A  2  2 

8(o*  >2<l-v  ) 


ttE 


f  \a  .  2  2,1/2  -  , 

__  ]  (a  -x  )  dx  J 


2 

J ~  (aa)  da 


(70) 


This  expression  yields  the  correct  result: 


2  A  2  2 

W  =  (oVJ  (1-v  )  /  E 
z  z 


(71) 


We  note,  however,  that  the  result  obtained  is  independent  of  the 
proportionately  factor  K. 

Now  suppose  we  try  the  approximate  trial  function 

AUz  =  K  (a  -  | x | )  | x |  <  a  (72) 


4U2  =  0 


x  >  a 


In  this  case  Eq.  (67)  yields  the  result 


W  =  2 • 67a2 (oA  )2(1-V2)  /  E 

2Z 


or 


Wfc  .  .  /  W 
trial  correct 


0.85 


This  is  a  fairly  crude  result,  however,  the  trial  function  is  also 
a  fairly  crude  one  relative  to  the  correct  choice.  But  with  this 
choice,  AU  has  dropped  to  half  its  maximum  value  at  x  =  0.5,  while 
the  true  function  has  dropped  to  half  its  maximum  value  at  x  =  0.87. 
So  the  stored  energy  should  be  higher  for  the  true  value  of  AU. 


CONCLUSION 

We  have  derived  a  variational  principle  to  determine  the 
value  of  the  stored  energy  or  acoustic  scattering  coefficient  of 
a  plane  crack.  We  use  the  crack  opening  displacement  as  a  trial 
function.  A  simple  example  for  a  plane  crack  using  a  triangular 
trial  function  for  AUz  gave  a  result  with  an  error  of  15%.  This 
would  lead  to  an  error  of  8%  in  estimating  K^.  The  use  of  Rayleigh 
Ritz  techniques  or  a  better  choice  of  trial  function  could  improve 
this  result. 

The  technique  is  very  easy  to  modify  to  take  account  of 

non-uniform  applied  stresses.  As  an  example,  if  we  take  a  bending 

stress  of  the  form  =  aB|x/a|,  AU  =  (a2-x2)^2,  we  find  that 

z  z 

I 

i 


W  =  0. 18TTa2  (oB)  2  (1-v2)  /  E 


This  would  imply  that  the  stress  intensity  factor  is  reduced  by 

0.42  from  the  value  with  a  uniform  applied  stress  aB.  The  result 

9 

for  the  exact  solution  is  a  reduction  0.5,  so  the  error  is  16%. 

Again  the  use  of  Rayleigh  Ritz  methods  would  be  helpful 
in  obtaining  more  exact  results. 
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